Positivity bounds on generalized parton distributions in impact parameter 

representation 



P.V. Pobylitsa 

Institute for Theoretical Physics II, Ruhr University Bochum, D-44780 Bochum, Germany 
and Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, 188350, Russia 

New positivity bounds are derived for generalized (off-forward) parton distributions using the 
impact parameter representation. These inequalities are stable under the evolution to higher nor- 
malization points. The full set of inequalities is infinite. Several particular cases are considered 
explicitly. 
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PACS numbers: 12.38.Lg 

I. GENERAL FORM OF THE POSITIVITY 
BOUNDS ON GPDS 

Generalized parton distributions (GPDs) [| [| g, |, 
^, fj], D also known as off-forward, skewed, nondiagonal 
etc. appear in the QCD description of various hard pro- 
cesses e.g. deeply virtual Compton scattering and hard 
exclusive meson production. Our knowledge about GPDs 
is poor and any additional theoretical information is of 
value. From this point of view the positivity bounds are 
rather important. 

GPDs arc defined in terms of matrix elements 



(x 1 ,X 2 ) = (P2,i2\4>L (^2) <j>cn (Xl) |Pl,«l> 



XoiCM^^Olo). (1) 



-{P2,i 2 \Pi,h 



Here \Pk,ik) is a hadron state (with momentum Pk and 
spin/isospin indices ik) and the field <p ak (xk) describes 
the annihilation of a parton with momentum fraction Xk 
and with spin/isospin labeled by a^. Indices also 
contain information about the type of the parton (quark 
or gluon). The vacuum subtraction term in the RHS of 
Eq. (Q) can be ignored in practical applications of GPDs 
with P\ 7^ P 2 but this term is important in the derivation 
of the positivity constraints on GPDs. Parton momen- 
tum fractions Xi = /P + are normalized with respect 
to some external fixed scale P + and not to the hadron 
momenta Pi or P 2 . Although notation ([!]) for GPDs dif- 
fers from the standard one, we find the form ([!]) rather 
convenient for the derivation of positivity bounds and for 
the analysis of the interplay between the evolution and 
the positivity properties. 

The positivity of the norm in the Hilbert space of states 



f dP + d 2 P ± dx , s , . , s , ( 
E J 2P+(27r)3 6(x)f ia (x,P)ct> a (x)\P,i) 



> 



leads to the inequality 

^ J dP+d 2 P^dxi f dP+d 2 P 2 ± dx 2 

i 1 a 1 i 2 a 2 



(2) 



2P+(27r)3 



2P+(27T)3 



0{xi)6(x 2 ) 



x f* 2a2 (x 2 ,P 2 )f lia Axi,Pi)G a P lt:, P2l2 (xi,x 2 ) > 0. (3) 
The integration over 2:1,2 is restricted to the positive re- 



gion 



X\ > , x 2 > 



(4) 



where the vacuum term vanishes in the RHS of Eq. (|l|) 
so that its subtraction does not violate the positivity. In 
the case of the antiquark GPDs one should consider the 
region xi,x 2 < 0. 

Strictly speaking the above expressions have to be writ- 
ten more carefully: in the case of the gauge theories the 
P exponent should be inserted between parton fields and 
certain conventions have to be chosen concerning the nor- 
malization of the momentum fractions xi,x 2 . 

Various positivity inequalities for GPDs corresponding 
to specific choices of the function fiajx, P) have been al- 
ready discussed in literature §, |l(J y, g |l|, [b| fTJ| 
[Tt], [l8| . The aim of this paper is to analyze inequality 
(0) with arbitrary functions fi a (x, P). 



II. POSITIVITY BOUNDS AND EVOLUTION 

It is well known that in the case of forward parton 
distributions the probabilistic interpretation of the one- 
loop DGLAP evolution Q |2l], |§ leads to the stability 
of the positivity properties under the one-loop evolution 
to higher normalization points. The generalization of this 
property for some particular positivity bounds on GPDs 
was considered in Ref. [JO]. Let us show modifying the 
argument of Ref. |Tl[| that if at some normalization point 
inequality (|3|) holds for all functions fi a (x,P) then after 
the one-loop evolution 



d_ 
dfx' 



a— r< a i a 2 . ,. 

a* p > . u p 1 i 1 ,p 2 i2 y xi 



(xi,x 2 ;n) = g 2 {y) / dyidy 2 



Pi 02 



ft£fl.i 3 (»i'»»5M) (5) 



to a higher normalization point inequality (g) is still valid 
for all f ia (x,P). 



Positivity bound 



involves 



GPDs 



G p 1 1 i 1 2 ,p 2 i 2 ( Xl > X2 >t i ) witn positive xi,x 2 > (Q). 
If x 1; X2 > then the one-loop evolution kernels 
Kp*0* (x\ ,X2,yi, 2/2) differ from zero only in the region 



Vk > x k > . 



(6) 



This constraint has a simple physical meaning: a par- 
ton with momentum fraction can emit a parton with 
momentum fraction Xk only under the condition (Q). 

The one- loop evolution kernels K^p* (xi, X2, yi, 2/2) 
for GPDs can be interpreted as perturbative parton-in- 
parton GPDs. For parton-in-parton GPDs one can re- 
peat the derivation of the inequality (^) arriving at the 
following inequality holding for any function h a p(x,y) 
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FIG. 1: Longitudinal momentum flow corresponding to vari- 
ables x,£ [in units of (P^ + P 2 + )/ 2 ]• 



Let us introduce a light-cone vector n and define the 
light-cone coordinates so that (nX) = X + for any 4- 
vector X . Below for simplicity we shall restrict our anal- 
ysis to the case of "unpolarized" GPDs defined in terms 
of the matrix elements of the operator ip(n'y)ip 



dX 
2^ 



exp 



iX(kf + k+) 



xh ollPl (x 1 ,yi)K^^(xi,x 2 ,yi,y 2 ) > 0. 



(7) 



The diagrammatic interpretation of this inequality is de- 
scribed in Appendix 

Actually inequality (R) holds only for functions 
hapix, y) vanishing at x — y because of the virtual terms 
proportional to 5{x\ — yi)6(x2 — 2/2) which give a neg- 
ative contribution to the kernel Kg^g'(xi, X2, 2/1, 2/2)- If 
one ignores these terms then the stability of the positiv- 
ity bound (^) with respect to the evolution (^|) upwards 
in {j, is a consequence of the property (]?]) of the one- 
loop evolution kernel. The last step is to notice that 
after the inclusion of the virtual terms proportional to 
<5(xi — yi)S(x 2 — 2/2) in the evolution equation, the posi- 
tivity is still preserved under the evolution upwards. One 
can use the same argument as in the case of the forward 
distributions: if at some "critical point" /x inequality (|^) 
is saturated and becomes an equality for some function 
fi a (x,P) then at this point the virtual terms propor- 
tional to 5(x% — yi)5{x 2 — 2/2) do not contribute to the 
evolution of the LHS of inequality (^) so that at higher 
\x the positivity is again restored by the positive part of 
the evolution kernel. The details can be found in Ap- 
pendix H. 
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(8) 

The structure of the momentum dependence in the RHS 
is fixed by the Lorentz invariance. In order to see the 
relation between the arguments of function P CT2Cr i and 
the standard notation of X. Ji |12| 



A = P 2 - Pi , t 



, 9 (nPi) - (nP 2 ) 

a 2 , £ = : (9) 



(nP 2 ) + (nPi) 



one can use the frame where P 2 



A -72. In this 

frame the RHS of Eq. (g) is simply F a2<7l (x, £, A 4 -). The 
longitudinal momentum flow corresponding to variables 
x, £ is shown in Fig. ^ 

Let us insert Eq. (g) into the positivity condition (||) 
and turn to the impact parameter representation 



Fa 2 o\ (x, b 



d 2 A A 
(2nf 



exp [i{AH x )] 



III. POSITIVITY BOUNDS IN THE IMPACT 
PARAMETER REPRESENTATION 

The positivity bound (|^) contains a multidimensional 
integral in the LHS and this form is not quite convenient 
for practical applications. Turning to the impact param- 
eter representation jl4| [l?], [l^, [l9| one can simplify this 
inequality. Let us show how this can be done for quark 
GPDs. The generalization for gluons is straightforward. 



xP CT2ffl (x,£,A x ) 



(10) 



Using a factorized ansatz for the function / ilQl (xi,Pi) 
one can reduce inequality (||) to the following rela- 
tively simple form (technical details can be found in Ap- 
pendix 



12 I I dx (i- x f p ^ (^-7) pi 1 1 



3 



> 0. 



(11) 



This inequality should be valid for any functions p a (orig- 
inating from the factorized ansatz for fi iai (x%, Pi)) and 
for any value of parameter b 1 - if one wants the original 
inequality (Q) to hold for arbitrary fi lUl Pi)- 

For hadrons with spin (e.g. pions) we define the GPD 
as follows 

dX 



q{x,£,t) = / — exp(iAx) 

Z7T 



■ \~) {ni)v (*£) Pi) 



(12) 



and for hadrons with spin 1 /2 we use the standard nota- 
tion of X. Ji || 

^exp(i\x){U(P 2 M (WW (t) |C/(Pl)> 

= H(x,£,t)U(P 2 )(n~f)U(Pi) 



IV. POSITIVITY BOUNDS ON NUCLEON 
GPDS 

From the practical point of view the most interest- 
ing case corresponds to the positivity bounds on nucleon 
GPDs. Let us rewrite our general inequality ([ll]) for 
spin-1/2 hadrons in terms of the standard notations for 
nucleon GPDs HJZ ©. 

Inserting Eqs. (10) and ( |l6|) into the general inequality 
(0) we find 



l€l 



1-^ 



ri 2 A J 



■ exp 



1-C : 



-(A^) 



2H 



2j* 



■P 



-A 1 +iA 2 j? 



A^iA 2 
A/(l-C 2 ) 



M(l-£ 2 ) 



E 2H- 



2j z 



P 



> 0. 



(17) 



+ 2M E(X ' C ' *) C/ ( p 2)^n M A l/ C/(P 1 ) 



(13) 



Here we assume the normalization condition n(Pi+P2) = 
2. Note that we could not impose this condition earlier 
since Pi and P2 were integration variables in Eq. (0) . We 
also follow the standard convention that the transverse 
coordinates are orthogonal to both n and Pi + P2 . The 
Lorentz invariant squared momentum transfer t = A 2 < 
can be expressed in terms of the transverse part lA -1- ! 2 



*=-(i-^r i (iA- 



| 2 +4£ 2 M 2 ) 



(14) 



In the case of spin-0 hadrons the function P CT2CTl ap- 
pearing in inequality ([ll]) can be expressed in terms of 
GPD q @ as follows 



F(x,£,A x ) = q [x,Z 



IA 



_L|2 



4£ 2 M 2 



(15) 



For spin-1/2 hadrons we have 



Pr 



2P- 



l (i,e,A- L ) = v/t^ 



2g- 



T P 



A^iA 2 
A/(l-? 2 ) 



-A'+iA 2 p 
M(l-£ 2 ) ^ 



P 2P - 



_2C 



tP 



(16) 



where the arguments x, £ and t (|1J) are assumed for H 
and P in the RHS. 

Inequality ( pd| ) is our main result. Combining the ex- 
pressions for Pcrjjo-! in terms of GPDs (|l5|), ( [l6|) with the 
impact parameter representation ( |l0| ) and inserting the 
result into inequality ( |ll| ) one can obtain the explicit 
form of the inequalities for GPDs q, H, P. 



The integration over A 1 - can be expressed in terms of the 
integration over variable t (|l4|). Since functions p a (z) are 
arbitrary we can rescale them: p a (z) — > z 3 / 2 p CT (z). Then 
inequality (|l7j ) takes the form 



1 1 



to 



(1-x) 2 



P+-P-,+ +P- 



dt Jq 



P-.-P-, 



p 



0> 



|6 x |M(l-a;) dt 
Here Jo is the Bessel function, parameter 

4£ 2 M 2 



[(t - t)E] > . (18) 



to 



1-e 2 



(19) 



corresponds to the maximal kincmatically allowed value 
of t (0) , functions P, P are taken at their standard ar- 
guments x, £, t, and p a ,± are the values of arbitrary func- 
tions Pa(z) at points z = (1 — ^ 



P<7,± = P<T 



1-X 



(20) 



V. SPECIAL CASES 

The general positivity bound ([ll]) imposes rather se- 
rious constraints on GPDs since it should hold for any 



4 



functions p„ and for any values of the impact parameter 
b . Let us show how choosing various functions p a one 
can reproduce most of the old positivity bounds and ob- 
tain new interesting results. For simplicity we consider 
the case of spinless hadrons (the generalization for spin- 
1/2 hadrons is straightforward) . 

1. Integrating inequality (O) over with the weight 
|£,-L|2n anc [ replacing p(z) — > z n+2 p(z) we obtain the fol- 
lowing inequality for GPD q ( |l5| ) 



l-x\ fl-x 



dt; J dx 

lei 



> 



t=-4£ 2 M 2 /(l-£ 2 ) 



(21) 



which should hold for any integer n = 0, 1, 2, . . . and for 
any function p. 

2. Taking p(z) = S(z — c) in (111]) one arrives at 



/ f^ eXp [ ?(A±6±) ] flfoMU-IA-T >0- (22) 

This coincides with the inequality derived earlier in Ref. 
@- 

3. Using function 



p (z) = aS (z - Y~^j + C 2 S - 
in inequality ( fTT| ) we obtain 



(23) 



|C2 



; j ■ v \ I X — £ 1 — X , 



i-r 0, i-r 



^ 2 + c^/ u-* &J 



> 0. 



(24) 



Taking into account that this inequality should hold for 
anyci,c 2 and rescaling b- 1 — > (1 — £ 2 )(1 — x)~ 1 b ± we find 



< JF(n,0, (1 - 0^)^(^2,0, (1 + 0b x ) 



where 



H, 2 = (a;±0/(1±0- 



(25) 



(26) 



This inequality was derived in Ref. [^7| (with a different 
normalization used in the definitions of GPD F and of 
the impact parameter b ). 

4. Taking Cfe = dk exp [—iP^b 1 -^ in inequality (p4]), in- 
tegrating over b and optimizing the resulting inequality 
with respect to arbitrary coefficients dk we obtain 



< (a;i)q(x 2 ) 



(27) 



where x±,X2 are given by Eq. ( J2q ) and #(0:1,2) = 
9(^1, 2> 0, 0) is the usual forward parton distribution. This 
inequality was obtained earlier in Ref. |ll| (the authors 
of |ll| use a different normalization of GPD q ). 

5. Using the modification of the ansatz ( |23j ) corre- 
sponding to spin-1 12 hadrons one can derive from the 
general inequality (pi]) 



H 



e 



i-e 



-E 



E 



2M(1 - 



< 



q{x 1 )q(x 2 ) 



(28) 



Again relation ( pq ) is assumed between the arguments 
x,£ in the LHS and the variables 2:1.2 in the RHS. This 
inequality was derived earlier in Ref. 

This short list of particular cases is only a small part 
of the bounds that can be extracted from the general 
inequality (|TT|) . 



VI. POSITIVITY BOUNDS AND 
RENORMALIZATION 

GPDs are defined in terms of matrix elements ( |1^ 
( |l3| ) of parton fields separated by a light-cone interval. 
Similarly to the case of the forward parton distributions 
(FPDs) this formal definition makes sense only in a com- 
bination with some renormalization procedure. Gener- 
ally speaking the renormalization includes subtractions 
which can violate naive positivity bounds. 

Working with the regularizations preserving the pos- 
itivity of the norm in the Hilbert space of states (with 
all necessary comments concerning the color singlet sec- 
tor, the insertion of the P exponent between parton fields 
etc.) one seems to be on the safe ground, which gives an 
argument in favor of the validity of the positivity bounds 
at high normalization points. On the other hand, in the 
case of FPDs it is well known that only the cross sec- 
tions associated with FPDs must be positive whereas 
the naive positivity bounds may be violated for FPDs 
at low normalization points in nonphysical renormaliza- 
tion schemes p^j. One should keep in mind that start- 
ing from the general inequality (|3j) one can reproduce 
the standard positivity properties of FPDs. Therefore 
a violation of the positivity properties of FPDs at low 
normalization points would lead to the breakdown of the 



5 



z r 



FIG. 2: Diagram representing the evolution kernel 
K^^(x 1 ,y 1 -x 2l y 2 ) at < x k < y k . 



positivity bounds on GPDs (this can be directly seen in 
inequalities (p7|), ( p8| ) where the GPDs are constrained 
by FPDs). 

In this paper we have shown that the validity of the 
positivity bounds on GPDs at high normalization points 
is compatible with the one-loop evolution. This self- 
consistency check is encouraging but certainly more se- 
rious analysis is needed in order to clarify the status of 
inequality (^) in the context of the renormalization. 



In the framework of Ref. |24| this contribution naturally 
appears from the pole integration over the light-like com- 
ponent of the momentum. In the case xi,X2 > the 
poles of the propagators associated with x\ and X2 lie 
on the same side of the real axis. Shifting the integra- 
tion contour to the opposite side one can get a nonzero 
contribution only from the pole of the propagator corre- 
sponding to the intermediate parton in the diagram of 
Fig. |2|. Functions V ak j3 kl {xk, yk) correspond to the two 
vertices of the diagram of Fig. The kinematical and 
normalization factors associated with the cut propaga- 
tor are obviously positive (this positivity is explicit in 
the light-cone gauge) and can be in clud ed into factors 
Va k f3 k -y(xk,yk)- In the RHS of Eq. (Al) one sums over 
the polarization 7 of the intermediate parton. Using the 
light-cone gauge and the techniques of Ref. it is easy 
to see that the summation over 7 is restricted to the phys- 
ical polarization s of the intermediate parton. Note that 
we can rewrite ( Al ) in the form 



^lift?^ 1 ' ^2; 2/1,2/2) 



VII. CONCLUSIONS 

Using the impact parameter representation we have 
derived positivity bound (|ll| ) on GPDs. This inequal- 
ity should hold for any function p,y so that actually we 
deal with an infinite set of inequalities. These positivity 
bounds impose certain constraints on models of GPDs 
used in the phenomenological analysis of hard exclusive 
processes. 
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APPENDIX A: POSITIVITY PROPERTIES OF 
EVOLUTION KERNELS 

In this Appendix we show how the positivity prop- 
erty ^ of the evolution kernel (xt, x 2 ;y 1,2/2) can 
be seen in the direct diagrammatic calculation of these 
kernels. Generally speaking, the leading order evolution 
kernel gets contributions from several diagrams |p3 but 
in the region < Xk < yk (||) which is interesting for us 
this kernel is given by the single cut diagram of Fig. ||. 
This diagram leads to the following structure of the evo- 
lution kernel 



(^1^2; 2/1, 2/2) 



E 

7 



<%i -xi-y 2 + x 2 ) 



x v ai ^(x!, yi)V* 2 p 2J (x 2 , 2/2) 



(Al) 



dz 
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VF Qlft7iZ (a;i,yi)VF* 2ft7 ^:E2,2/2) (A2) 



where 



W a p~f, z (x,v) =S(y -x- z)V af 3~ ( {x,y) . 



(A3) 



Obviously the form ( |A2| ) of the evolution kernel K au- 
tomatically leads to the positivity property (|t]) . Actu- 
ally the general decomposition ( |A2j ) is essentially equiv- 
alent to the inequality (0) . It is decomposition ( |A2[ ) that 
is needed for the proof of the stability of the positivity 
bounds on GPDs with respect to the evolution to higher 
normalization points. 

In order to illustrate the above general formulas with 
an explicit example let us consider the evolution kernel 
for the "helicity-independent" quark GPDs 



K 0il3 2 ( x i> x ^yi,y2) = S/3 lf 3 2 K(x 1 ,x 2 ;y 1 ,y 2 ) . (A4) 



In the region (||), < Xk < yk, the kernel K is given by 
the following expression p4[ 



K(x 1 ,x 2 ;y 1 ,y 2 ) 



2tt 2 



HVl ~ X l ~ 2/2 + X 2 ) I t - ' ' J - ' ' ;2 



2/1 - xi 



2/12/2 



(A5) 



Obviously this kernel has the form (Al) 



K{xi,x 2 \yx,V2) 



Cf_ 

2tt 2 



5(yt - xi - y 2 + x 2 ) 
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1 



1 



X 2 



vm - xi \Jv2 - x 2 yw/yi - %i vi\Jv2 - x 2 

(A6) 

This contribution comes from the diagram of Fig. m with 
a gluon playing the role of the intermediate parton. The 
two terms in the brackets in the RHS correspond to two 
possible polarizations of this gluon. 

The above expression for the kernel K is singular at the 
points yk — %k ■ The proper treatment of these singulari- 
ties is described by the following expression for the con- 
volution of the kernel K with arbitrary functions h(x, y) 

S 1 n 



dxi / dyi / dx 2 / dy 2 h(xi,yi)h*(xi,yi) 



OC OO i OO 

:if(:ri, £2; 2/1,2/2) = ■^lim J dxi J dx2 \ J ~ 



1 + 



£12:2 



(xi + z)(x 2 + z) 



h(xi,xi + z)h*(x 2 ,x 2 + z) 



3 \.r..r> , 

- - In— 2— ) h(xi,xi)h (x 2 ,x 2/ 



(A7) 



The last term in the RHS proportional to 
h{x\, Xi)h* (x 2 , x 2 ) corresponds to the virtual 
5{y\ — xi)6(y 2 — x 2 ) contributions to the evolution 
kernel K . In the limit e — > both integral and contact 
terms are divergent but these two divergences cancel 
each other. 



APPENDIX B: STABILITY OF POSITIVITY 
BOUNDS UNDER EVOLUTION 

In this Appendix we present a detailed derivation of the 
stability of the positivity bounds under the one-loop 
evolution (^) to higher normalization points. We shall 
use compact Dirac notation for the integral appearing in 
the LHS of inequality (^) 

if\r\n- V f dp i +d2p i ±dx i f dP 2 + d 2 P 2 A -dx 2 

tfM/;- Z_. J 2P+(27r)3 J 2P+(27r)3 



x0( Xl )e(x 2 ) /* a2 (x 2 ,P 2 )f liai (xx,Px)G%% P2i2 ( Xl ,x 2 ) 

(Bl) 

so that the positivity bound (|3j) can be written as follows 



Let us assume that at some normalization point in- 
equality (H) holds for all functions fi a (x,P). Imagine 
that during the evolution to higher /1 this inequality 
breaks down at some point /io for some function so 
that at \i > /io 

(/ (0) |G M l/ (0) )<0 (^>Mo) (B3) 
but at (i < /io we still have for all functions / 

(/|GJ/)>0 (/i<Mo). (B4) 
Then at point /j,q 

</ (O) |GJ/ (o) )=0. (B5) 

The fate of this "matrix element" at /i > /io is determined 
by the evolution equation (||) 

^(f {0) \G„\fW) = (fW\(K®G fl )\fM). (B6) 
a/i 

Here we use short notation K ® G M for the convolution of 
the evolution kernel K with GPD G. If one could show 
that at the point /io we have 



{f^\{K®G m )\f^)>Q 



(B7) 



then this would guarantee that the evolution from /to to 
higher /i would restore the positivity of (f^lG^lf^). 
This would invalidate our assumption (|B3j). Thus in or- 
der to prove the stability of the positivity bounds on 
GPD under the evolution to higher normalization points 
we need the following 

Statement. If at some normalization point /io for all 
functions / 



and for some Z^ -* 



(/|G M0 |/> >0 



(J (0) |GJ/^) =0 



then 



(/W|(^®Gpo)|/W)>0. 



(B8) 



(B9) 



(BIO) 



To prove this statement we make use of the positivity 
property of the evolution kernels (Q). As it is explained 
in Appendix ^ this po sitivity property is equivalent to 
the decomposition ( |A2j ) of the evolution kernel. In our 
short Dirac-like notations Eq. ( |A2| ) takes the following 
form 

K = j dzJ2\W 7 , z )(W 1>z \ (Bll) 



Inserting repr esentation ( Bll ) for K into the LHS of in- 
equality (BIO) we find 



(f\G\f) >0. 



(B2) 



(f^m^G^ifW) 
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dzJ2(fW (^W^lG^lfW ®W 1>Z ) >0. (B12) 



The RHS is positive according to Eq. (|Bq) . This com- 
pletes the derivation of the inequality ( |B10| ) . 

The above proof of the stability of the positivity 
bounds ([3]) under the evolution (|5|) ignored the problem 
of the negative terms proportional to 8(x\ — yi)5(x 2 — y 2 ) 
which are present in the evolution kernel K. Because of 
these terms it is allowed to use inequality ([?]) only for 
functions hjx. y) w hich vanish at x = y. As a result in- 
stead of Eq. ( |B 1 1 ) one has to work with the following 
representation for the evolution kernel 



K 



lim (k^ 

ho 



id 2 ) 



(B13) 



where the first term is given by the e regularized integral 
decomposition ( Bll ) 



X2,yi,V2) 



*(xi,yi)W* tZ (x 2 ,y 2 ) 
(B14) 

and the second piece 

KP(xi,x 2 ,yi,y2) = C(e,xx,x 2 )S(x 1 - yi)S(x 2 - y 2 ) 

. (B15) 

comes from the virtual contribution to the evolution ker- 
nel which is proportional to 8{x\ — y\)5{x 2 — y 2 ) with a 
divergent coefficient C(s 7 Xi,x 2 ). This coefficient is reg- 
ularized by the same small parameter e so that 

C(e,xi,x 2 ) = cx + c 2 \ne + F 1 {x 1 )+F 2 (x 2 ). (B16) 

The singular contribution lne compensates the diver- 
gence of the z integral in th e RHS of Eq. (B14). Ob- 
viously representation 



513| ) guarantees that K obeys 
inequality (^) for arbitrary functions h(x,y) vanishing at 
x = y. An example of the gener al st ructure ( B16| ) can be 
seen in the explicit expression (A7) for the quark evolu- 
tion kernel where the term ln((4xix 2 ) / e 2 ) represents all 
terms appearing in the rhs of (B16). 



Now we hav e to s how that the kernel K (B13) still has 
the p rope rty ( |B1C| ). For the first term in the RHS of 
Eq. (B13) this can be done in the same way as in the 
inequality ( B12| ): 



(/ (0 W> ® G^l/W) >0. 
Next let us show that 

(f^\(Kp®G, )\f^)=Q. 
According to Eqs. ( pl5| ), ( |B10| ) we have 

(f^\(K^ ® G M0 )|/W) = </<°>|[( Cl +c 2 ln £ )G, 



(B17) 



(B18) 



Note that \f^) is associated with /W ( xi, y i) whereas 
(/(o)| corresponds to f (0) {x 2 ,y 2 ) in Eq. @. This ex- 
plains the ordering of "operators" F\ and F 2 associated 
with fu nction s F\(xi) and F 2 (x 2 ) respectively in the RHS 
of Eq. (B19). Taking into account Eq. (B9) we get rid 
of the ci + c 2 In e contribution: 



-(2) 



(0)\ 



</^|(G WJ F 1 + F 2 G /Jo )|/ 



(o)\ 



(B20) 



Since property of G W) (B8) holds for any / we conclude 
from Eq. (H) that 



G w |/W) = 0, (/(°)|G M =0. 



This means that 



(f ( - 0) \(G w F 1 + F 2 G lio )\f 



(0) 



(B21) 



(B22) 



Inserting this result into Eq. ( B20| ) we obtain Eq. (B18). 
Combining Eqs. pT?| ), ( |B18| ) we derive inequality ( |B10[ ) , 
which means that the virtual terms in the evolution ker- 
nel do not violate the stability of the positivity bounds on 
parton distributions under the evolution to higher nor- 
malization points. 



APPENDIX C: DERIVATION OF POSITIVITY 
BOUNDS IN IMPACT PARAMETER 
REPRESENTATION 

In this Appendix we show how inequality (^|) contain- 
ing a multidimensional integral can be reduced to the 
relatively simple form (pT[). 

Using representation (g) for GPD G we find from in- 
equality 



00 00 00 

Jdk+J^j-J <fpt J dJ4 



p+ 



d 2 pi 



2 -kt)-{Pt -kt)] 



k+ + k+ p+-p+ 2(p+p£-p+pj-y 



■}+ ' 



1 2 1 1 



Pi 



Po 4 



pi 



pi 



xfZ 2 {k+,P+,P.i)f ai (k+,Pt, Pi) > . 
Taking functions f a in the factorized form 



(CI) 



+G l _ l0 F 1 + F 2 G fJi0 ]\f 



(0)\ 



(B19) 



/, t( /, + .P~.P- l .=</„(A + .P+), / (| t 



(C2) 



8 



we derive from inequality (CI 



oo oo oo oo 



J dk+ J dp+ J dk+ J dp+d* 2 (k+,p+)d al (k+,p+) 







In the limit 

\g(u)\ 2 -> 6{u-u) 
we find from inequality (|C3| ) 



(C5) 



k-^ ~\~ k<2 



P l- P 2 P l+ P 2 



K+ P 2 + P l+ P 2 ^ P 1 P 2 



'S[(P+-k+)-(P+-kt)] 



P t+ P 2 + 
2P+P+ 



>0 (C3) 



where F a2<yi is defined by Eq. Jl(i|). This inequality 
should hold for any value of b 1 - and for any function 
d a (k + , P + ). The next step is to take 



d a {k + ,P+) = g{P+ - k+)h a (P^ 



(C4) 



l l 
d£ J dx 

- 1 iei 



n /',-,-, ( F a2ffl ( Xl e, %-4St ) > • ( C6 ) 



1 - X 



Ki-^ 2 ) 



Since function ha and parameter 6^ are arbitrary we can 
set v — 1 without losing generality. Taking h a (z) — 
z 2 p a -(z~ 1 ) we obtain inequality (11). 
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